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I. INTRODUCTION 

Based on earlier ideas of Wigner ^ Dyson introduced 
a three-fold classification of quantum systems accord- 
ing to their behavior under time-reversal and spin ro- 
tation 0. This symmetry classification turned out to 
be very useful, for instance in semiclassical, disordered 
and random-matrix approaches to complex quantum sys- 
tems. The success of random-matrix theory is based on 
universal features in spectra of complex quantum sys- 
tems. While not capable of predicting single eigenval- 
ues random-matrix theory has become one of the key 
ingredients in predicting physical features that depend 
on non-trivial spectral statistics Q . In each symme- 
try class various universality classes have been identi- 
fied - each described by some ensemble of random ma- 
trices. Most prominent are the three Gaussian ensembles 
QUE, GOE and GSE. They define the ergodic universal- 
ity classes and they have been applied successfully to a 
wide range of quantum systems (see the recent review jlj 
for an overview and further references). 

Recently the three-fold classification has been ex- 
tended to a ten-fold classification by incorporating spec- 
tral mirror symmetries. This lead to seven novel 
symmetry classes 0, 0, 0- They are partly real- 
ized for a Dirac particle in a random gauge field, 
and for quasi-particles in disordered superconductors or 
normalconducting-superconducting hybrid systems. In 
the presence of a spectral mirror symmetry the spectrum 
is symmetric with respect to one point Eq: ii Eq + E is 
in the spectrum so is Eq — E. The invention of the novel 
classes has become necessary due to the impact such a 
symmetry has on spectral correlations. These new uni- 
versal features appear near the symmetry point Eq and 
they can be described by defining random-matrix ensem- 
bles which incorporate the corresponding spectral mirror 
symmetry. 
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It has been conjectured by Bohigas, Giannoni and 
Schmit that the spectra of classically chaotic systems 
display the spectral fluctuations described by the three 
Gaussian Wigner-Dyson ensembles of random-matrix 
theory Though the fidelity to the universal pre- 

dictions of random-matrix theory have an overwhelming 
support by both experimental and numerical data the 
physical basis of universality is not completely under- 
stood. 

Quantum graphs have been introduced by Kottos and 
Smilansky ,1(1] as simple quantum models with an ex- 
act semiclassical trace formula for the density of states 
which is expressed as a sum over periodic orbits on the 
graph. They have since become an important tool in the 
semiclassical approach to universality. In this series of 
papers we will construct star graphs for all ten symme- 
try classes and investigate their spectral statistics both 
numerically and analytically. While the following paper 
[Tl| is devoted to a semiclassical periodic-orbit approach 
this paper focusses on the construction of appropriate 
star graphs and some numerical results. 

We start with giving a short introduction to the ten 
symmetry classes in section with all details needed 
for the subsequent construction of star graphs. In the 
following section IIIII on spectral statistics we introduce 
the spectral form factors, review the results of Gaussian 
random-matrix theory for the ten symmetry classes and 
generalize the Bohigas-Giannoni-Schmit conjecture. Af- 
ter a general introduction to quantum graphs in section 
II VI we construct one ensemble of star graphs for each of 
the ten symmetry classes. Numerical results then show 
the fidelity of these ensembles to the predictions of the 
Gaussian random-matrix ensembles. 



II. THE TEN SYMMETRY CLASSES OF 
QUANTUM SYSTEMS 

In quantum mechanics most symmetries are described 
by some unitary operators U that commute with the 
Hamilton operator Tl = UTiU^ . Thus the operators U 
(or its hermitian generators) describe constants of motion 
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and they lead to a block diagonal form of the Hamilton 
matrix in an eigenbasis ofU. If enough constants of mo- 
tion Ui are available such that the corresponding hermi- 
tian generators form a complete set of commuting observ- 
ables the Hamilton operator is eventually diagonalized 
in the common eigenbasis of the symmetry operators (or 
their generators). However, for any hermitian Hamilton 
operator Ti. there is always a complete set of commut- 
ing hermitian operators Vi which also commute with the 
Hamilton operator (e.g. projectors on eigenstates). In 
some sense the notion of symmetry in this wide sense is 
obsolete for a single quantum system. However, in most 
cases such a set of commuting operators will not have any 
corresponding classical observable and will only apply to 
a single system. 

It is more appropriate for our purposes to consider 
a family or class of quantum systems. Such a family 
may arrise by varying some physical parameters (like the 
strength of an applied magnetic field) or, for disordered 
systems, by an ensemble of random potentials. In the 
derivation and application of semiclassical methods one 
formally considers the asymptotics /i — > which is equiv- 
alent to a family of operators with fixed h but some vary- 
ing physical parameters. 

A unitary operator Ui is a unitary symmetry of a class 
of systems if it commutes with all Hamilton operators 
in that class. This notion avoids "symmetry" operators 
that would rely on a diagonalization of the Hamilton op- 
erator. In semiclassics a unitary symmetry will have a 
classical correspondence. In the sequel we will consider 
each Hamilton operator as a member of a class without 
explicitly referring to it. 

If a family of systems has a unitary symmetry all its 
Hamilton operators can be brought to a blockdiagonal 
form. Each block can be regarded as a new Hamilton op- 
erator on a reduced Hilbert space. Let us assume that the 
Hilbert space is completely reduced such that there are 
no more unitary symmetries. What types of symmetry 
may such a reduced quantum system still have? What 
are the possible structures of the Hamilton operator (or 
the Hamilton matrix) and what are the consequences on 
its spectrum and its eigenvectors? Such questions were 
for the first time addressed and partially answered by 
Wigner and Dyson P, Ql- Dyson proposed a symme- 
try classification based on the behavior of quantum sys- 
tems under time-reversal and spin rotation. This lead to 
three symmetry classes (the three-fold way): i. systems 
that are not time-reversal invariant, ii. time-reversal in- 
variant spin-less particles, and Hi. time-reversal invari- 
ant particles with spin s = i. Time- reversal symmetry 
has immediate consequences on the form of the Hamil- 
ton operator: spin-less particles can be described by real 
symmetric Hamilton matrices in a time-reversal invariant 
basis, while systems without time-reversal invariance do 
not have any canonical basis and the Hamilton matrix 
remains complex. The influence of the symmetry class 
on spectral properties such as level repulsion has been 
investigated extensively within the field of random ma- 



trix theory HHIii. We will give more details on random 
matrix theory in section ITTTl 

Recently the Wigner-Dyson symmetry classification 
has been extended to a ten-fold way by including all 
different types of mirror symmetries in the spectrum 
0, S I3- the presence of a spectral mirror symme- 
try every eigenvalue Eq -\- E has a partner eigenvalue 
a.t Eq — E (we will set i?o = in the sequel without 
loss of generality). Below we will describe the various 
ways a mirror symmetry may arise and be combined with 
time-reversal invariance. This leads to the seven novel 
symmetry classes. As shown in there is a one-to-one 
correspondence between Cartan's ten-fold classification 
of Riemannian symmetric spaces and the ten symmetry 
classes of quantum systems. We will use the conven- 
tion to adopt the names given by Cartan to the different 
classes of symmetric spaces for the according symmetry 
classes. 

The novel symmetry classes are partly realized for 
Dirac fermions in a random potential (the chiral classes) 
7] and partly for quasi-particles in mesoscopic super- 
conductors or superconducting-normalconducting (SN) 
hybrid systems (the Andreev classes). It is possible to 
construct much more general systems in the appropriate 
symmetry classes, e.g. two coupled spins or a general- 
ized version of the Pauli equation (which includes the 
Bogoliubov-de-Gennes equation as a special case), and 
quantum graphs. Quantum maps which incorporate the 
corresponding symmetries have been discussed recently 
[l^ . Due to their simplicity graphs will be the focus of 
this work. The following discussion of symmetry classes 
is summarized in table 



A. Time-reversal invariance 

Quantum systems obey generalized time-reversal sym- 
metry if there is an anti-unitary operator T - the gener- 
alized time reversal operator - that first, commutes with 
the Hamilton operator 



[n,T] = o 



and second obeys 



= ±1. 



(1) 



(2) 



Anti-unitarity implies i. anti- linearity T(a|f) + /3|i^)) = 
a*\TO + P*]^!^) and ii. {T£,\Tu) = . 

For time-reversal invariant systems T changes the di- 
rection of time when applied to the Schrodinger equation. 
Equivalently, when T is applied to the retarded Green's 
operator 



Q+{E) 



1 



E + ie-n 



(3) 



one gets 

Tg+{E)T-^ 



E-ie-H 



g+{E)^ ^g_{E), (4) 
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which is just the advanced Green's operator. 

Time-reversal symmetry also effects other dynamic op- 
erators - such as the unitary time evolution operator 



U[t) 



JHtjh 



(5) 



Scattering problems can often be described by some uni- 
tary operator iS(£') that connects incoming and outgoing 
states of energy E. Time-reversal invariance leads to 



rs{E)r-^ 



-U{-t) 
--S{E)\ 



(6) 
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These equations also define time-reversal symmetry for 
quantum maps. The transformation of the time devel- 
opment operator follows immediately from the condition 
on the Hamiltonian. In scattering problems S{E) can 
be related to a unitary combination of Green's functions 
- for definiteness consider S{E) = Q+{E)Q^{E)~^ — 
\ — 2t\Q+ (E) and equation © follows from the trans- 
formation ^ of G+{E). 

We have used the term generalized time-reversal oper- 
ator because T need not be the well-known conventional 
time-reversal operator. For a particle in IR'^ the anti- 
unitary conventional time-reversal operator obeys 
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where s is the particle spin. This conventional time- 
reversal operator obeys T^onv = 1 if the spin quantum 
number is integer s = 0, 1, 2, . . . , and %?onv — ~1 if 
spin is half- integer s = ^ , | , . . . . Thus the most relevant 
and simplest realizations are for spin- less (^onv — 1) ^^'^ 
spin-i (T^onv = ~1) particles. 

When a given quantum system is studied one should 
be aware that a generalized time-reversal operator may 
still exist which commutes with the Hamiltonian while 
the conventional time-reversal operator may not com- 
mute with Ti.. 

The consideration of time-reversal symmetries leads 
to three symmetry classes: either a system is not time- 
reversal invariant, or it is time-reversal invariant - in the 
latter case the time-reversal operator either obeys = 1 
or = —1. These classes have been called Wigner- 
Dyson classes and their impact on the form of Hamilton 
matrices and universal spectral features will be discussed 
further in section Hi C ll and lllll Additional spectral mir- 
ror symmetries lead to the novel symmetry classes to be 
discussed below. 

Kramers' degeneracy occurs in time-reversal invariant 
quantum systems with = — 1. If |x) is an eigenvector 
with eigenvalue E, then due to time reversal invariance 
\Tx) = T\x) is an eigenvector with the same eigenvalue 
E. It is straight forward to show that \Tx) is orthogonal 
to Ix) using the properties of the time reversal symme- 
try operator. This degeneracy is well known for spin-^ 
particles with conventional time-reversal symmetry. 



B. Spectral mirror symmetries 

A quantum system has a spectral mirror symmetry if 
the spectrum is symmetric: for every eigenvalue E > 
there is another eigenvalue —E < 0. In general, there 
may be some vanishing eigenvalues E = 0. We will now 
discuss the symmetry operators related to such spectral 
mirror symmetries. 

According to a theorem by Wigner any symmetry op- 
eration on Hilbert space is either represented by a uni- 
tary operator V or an anti- unitary operator C. Now take 
any eigenstate {v) such that Ti-lv) — Elv) - it is obvi- 
ous that spectral mirror symmetry demands that either 
nVlv) = -EVlv) or nCliy) = -EC\iy). This condition 
on an eigenstate leads eventually to the condition that 
the Hamilton operator anti-commutes with either a uni- 
tary or an anti-unitary symmetry operator 



[P,H].=0 



or 



[C,H], = 0. 



(8) 



As an additional condition any state should be invariant 
upto a phase factor when acted upon twice with V oi C. 
It can be shown that it suffices to consider 



-P^ = ±1 



±1. 



(9) 



If a system is not time-reversal invariant one can always 
choose = 1 while in presence of time reversal symme- 
try either ^ 1 or = -1. 

Spectral mirror symmetries relate Green's operators at 
energy E and —E 



pg+{E)r-' = -G-i-E) 
cg+{E)c-' = -g+{-E). 

For scattering problems this leads to 

VS{E)V-'^ =S{-E)^ 
CS{E)C-^ =S{-E) 

and for the time development operator to 

VU{t)V^^ = U{-t) ^U{t)^ 



(10) 



(11) 



(12) 



The seven novel symmetry classes are obtained by all 
possible combinations of a spectral mirror symmetry with 
time-reversal symmetry (with the additional requirement 
that \P,T] = or [C,T] = if both symmetries are 
supposed to hold). First, there are three novel symmetry 
classes that are not time reversal invariant: either there is 
a unitary operator with = 1 or an anti-unitary with 
= ±1. In time-reversal invariant systems one has 
both unitary and anti-unitary spectral mirror symmetry 
operators: if a unitary operator V anti-commutes with 
the Hamilton operator so does the anti-unitary operator 
C = VT. As = ±1 and = ±1 this leads to four 
symmetry classes that combine time-reversal symmetry 
with spectral mirror symmetry: if = 1 either = 1 
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-1 
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-1 


Sp{N) 


CI 


+1 


-1 


-1 


Sp{N)/U{N) 


BD (D) 








+1 


SO{N) 


Dill 


-1 


-1 


+1 


SO{2N)/U{N) 



TABLE I: The ten symmetry classes of quantum systems. If 
a symmetry class obeys time-reversal symmetry or a spectral 
mirror symmetry the entry ±1 in the corresponding column 
indicates if the symmetry operator squares to ±1. The entry 
indicates that the corresponding symmetry is broken. The 
last column gives the corresponding Riemannian symmetric 
space (of compact type). 



(p2 ^ 1) or = -1 {V^ = -1), if = -1 either 
C2 ^ -1 (p2 ^ 1) or C^ = l {V^ = -1). 

For historic reasons these seven classes have been split 
into two groups, the first group is given by the three 
chiral classes ~ the ones that have a unitary mirror sym- 
metry with = 1. Their importance has first been 
observed in investigations of Dirac fermions in a random 
gauge field where the spectral symmetry is related to 
chirality. For this reason we will call V a chiral symme- 
try operator though in general V need not be related to 
chirality. The four remaining classes have mainly been 
discussed in connection to mesoscopic disordered super- 
conductors or superconducting-normalconducting hybrid 
systems where the anti-unitary mirror symmetry is con- 
nected to electron-hole conjugation. For this reason we 
call C a charge conjugation symmetry operator, though 
again, in general C need not be related to charge conjuga- 
tion at all. Since Andreev reflection is a main ingredient 
in the dynamics of superconducting-normalconducting 
hybrid systems we will call these classes Andreev classes. 
The detailed discussion of these symmetry classes and 
their impact on universal spectral features will be dis- 
cussed in sections III C 21 III C 31 and IIIII 



C. Explicit form of scattering matrices for each 
symmetry class 

Time-reversal and spectral mirror symmetries restrict 
the form of Hamilton and scattering matrices due to the 
relations (Q), ®, © and ((TT|l . By choosing an appro- 
priate Hilbert space basis for each symmetry class the 
symmetry operators are represented by a simple ma- 
trix (combined with the complex conjugation operator 



for anti- unitary operators). These determine the explicit 
form of scattering matrices for each symmetry class. 

Note, that the following derivation of the scattering 
matrices depends on the choice of the basis. There are 
many choices for the Hilbert space basis in which the 
symmetry operators have a simple form. As a conse- 
quence many of the following identities are only valid 
in that special basis. Especially the "complex conjuga- 
tion operator" /C is defined with respect to a given basis. 
However, one may show that the bases chosen here can 
always be constructed from the general properties of the 
time-reversal and spectral mirror symmetries. Our choice 
of basis is biased by their later application to star graphs 
in section ITVI 

In addition, some symmetry classes have a further di- 
vision into subclasses. Though we will mention all sub- 
classes we will only give the scattering matrix in one of 
the subclasses. 



1. The Wigner-Dyson classes 

Quantum systems without spectral mirror symmetries 
belong to one of the three Wigner-Dyson classes A, Al 
or All. 

Class A contains quantum systems that are not time- 
reversal invariant. There is no preferred basis in Hilbert 
space and the scattering matrix S{E) may be any unitary 
N X N matrix. 

A time-reversal invariant system belongs either to class 
Al if = 1 or to class ^11 if = -1. 

In class Al there are time-reversal invariant bases such 
that T\i) = \i) for any basis state. In any such basis the 
time-reversal symmetry operator is represented by the 
complex conjugation operator 

Al: r = /C, (13) 

where the complex conjugation operator acts on a gen- 
eral state by complex conjugation of the coefficients 

'^E»=ia«N> =''CE^ia»l*> =E^i<l*)- ^ 

The condition (jSJ implies that a scattering matrix is 

represented by a unitary symmetric N x N matrix 

Al : S{E)=S{Ey. (14) 

For class All there is no time- reversal invariant basis. 
Instead, there are always bases in which the time-reversal 
symmetry operator is be represented by 

All: T^jc(^^ (15) 

where 1 is the N x N identity matrix. Hilbert space has 
even dimension due to Kramers' degeneracy. In such a 
basis the scattering matrix has the form 

All : SiE) = l'''^^'' ""'^^^ I (16) 
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with complex N x N matrices Xi that satisfy X2{E) = 
-X2{EY and X:i{E) = -X^^iEY and are further re- 
stricted by unitarity of S{E). 



2. The chiral classes 

A system with a spectral mirror symmetry connected 
to a unitary chiral symmetry operator V (with = 1) 
falls into one of the three chiral symmetry classes yllll, 
BDl or CII. 

Since — 1 the eigenvalues are either +1 or —1. In 
general, there will be p positive and q negative eigenval- 
ues. The number v — \p — q\ = 0,1,2,... distinguishes 
between different subclasses in each of the chiral classes 
(y is always even for class CII) . The integer v has impact 
on both the form of Hamilton (or scattering) matrices 
and on the spectral statistics. Because V relates states 
with positive energy to states with negative energy there 
are u vanishing energy eigenvalues due to the chiral sym- 
metry. 

We will focus on the subclasses with v = and set 
p — q = N ill classes AIII and BDl, p = q ^ 2N in class 
CII. Hilbert space has even dimension in all three classes. 
There are many bases that can be used as reference basis, 
for example the one, where V is diagonal. Here, biased 
by our following construction of star graphs we choose 



AIII, BDl, CII 



V 




(17) 



which can be obtained from the diagonal representation 
by a simple rotation. 

The chiral class yllll contains systems without addi- 
tional time-reversal invariance. The other two chiral sym- 
metry classes are time-reversal invariant with ~ 1 for 
class BDl and T'^ — —1 for class CII. In class BDl one 
may always choose a time-reversal symmetry operator of 
the form 



BDl 



T = 1C 



(18) 



which commutes with the chiral symmetry operator V 
(I17|l . In class CII one may choose 



CII 



T = /C 



/o 


-1 





o\ 


1 




















-1 







1 





(19) 



which also commutes V. 

Due to the condition 1)11(1 a scattering matrix S{E) in 
class j4III has the form 



Alll : S{E) = 



Xt{E) X2{E) 
,X3{E) Xi(~E)'' 



(20) 



where (besides unitarity) the N x N matrices Xi are 
further restricted by X2{E) = X2{-Ey and X^iE) = 
Xsi-E)^. 



In class BDl, due to time-reversal invariance JHl, S{E) 
is symmetric, thus 



BDl : S{E) 



(21) 



CH: S{E) = 



(22) 



Xi{E) X2{E) \ 
^X2iE)'^ Xi{~Erj 

where Xi{E) = Xi{E)^ and X2{E) = X2{-Ey. 

Finally, in class CII, S{E) is a 4A^ x AN matrix of the 
form 

/ Xi{E) X2{E) XsiE) X4iE)\ 
Xs{E) Xi(E}'^ Xi{-E)^- Xe(E) 
XeiE)'^ -^XiiE)"^ Xii-E^ Xr^i-E^ 
\~X4{-Er X3(E)'^ X2(-E)'< Xi(^E)j 

with additional constraints X2{E) = —X2{EY', X^{E) = 
Xs{-EY, X^{E) = ~X.,{EY and X^{E) = X^-E)^ . 

3. The Andreev classes 

A quantum system with a spectral mirror symmetry 
that does not belong to any of the chiral symmetry classes 
belongs to one of the four Andreev classes C, CI, BD or 
Dill. The spectral mirror symmetry for these classes is 
related to an anti-unitary charge conjugation operator 
C with = -1 for C and CI while = 1 for BD 
and DHL The classes C and BD are not time-reversal 
invariant while CI and I?III are time-reversal invariant 
with = 1 in CI and = -1 in Dill. 

The classes C and CI do not split into subclasses. In 
appropriate 2A^-dimensional bases the charge conjuga- 
tion operator can be represented as 



C, CI : 



C = IC 







(23) 



In contrast the classes BD and Dili fall into two sub- 
classes each. The symmetry class BD allows for either 
an even or odd dimensional Hilbert space. Due to spec- 
tral mirror symmetry there is always an eigenvalue on 
the symmetry point £' = in an odd- dimensional Hilbert 
space. The subclass with odd (even) dimensional Hilbert 
space may be called i?£'-odd(even). In the following we 
will restrict ourselves to the even-dimensional case and 
will follow the convention to call it symmetry class D. 
Similarly Dill falls into the two subclasses I?III-odd and 
Dlll-even. The dimension of the corresponding Hilbert 
spaces is twice an odd or twice an even number. Spectral 
mirror symmetry combined with Kramers' degeneracy 
implies two eigenvalues = on the spectral symme- 
try point in class DIITodd. In the sequel we will restrict 
to Dlll-even which is physically more relevant. 

An appropriate choice of basis in the Hilbert space 
takes the charge conjugation operator of the symmetry 
classes D and Dili (we will not mention the "even" fur- 
ther) to the form 



i:>,L>III : C 




(24) 
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where 1 is the N x N identity matrix for class D and the 
2N X 2N identity for class Dili. 

The time-reversal symmetry operators in the classes 
CI and Dili have the representations 



CI: r = /c 



(25) 



and 



Dili : r = /C 



/O -1 o\ 

10 



Vo 1 / 



(26) 



in an appropriate basis - here 1 is the NxN identity ma- 
trix in both equations. These representations commute 
with the corresponding representations of the charge con- 
jugation operators. 

The conditions © and Hll|) to scattering matrices 
S{E) of the form 



C, CI ; SiE) 



Xi(E) X2{E) 
-X^i-EY Xii-E)' 



(27) 



in classes C and CI. There are no further restrictions 
on the complex NxN matrices Xi for class C (apart 
from unitarity). Time-reversal invariance in class CI re- 
quires S{E) to be symmetric, thus Xi{E) — Xi{E)'^ and 
X2{E) = -X2{-EY. 

In the symmetry class D the scattering matrix has the 
form 



D : S{E) = 



Xi{E) X2{E) 
,X2{~EY Xi{~E)\ 



(28) 



without further restrictions on the NxN matrices Xi. 

For class _DIII S{E) is a complex AN x 4N matrix of 
the form 



Dili: S{E)^ 



/ Xi{E} X2(E) X3(E) Xi(E)\ 

XrJE) Xi(Ef Xe{E) Xsi-E^ 
Xsi-E)* Xii-E)* Xii-EY X2{-Ey 
^ei-E)' Xsi-E)"^ Xsi-E)' Xii-EYj 



(29) 



with X2{E) = -~X2{E)'^, Xi{E) = 
-X^{Ef and X7{E) = -X-j{-E)K 



^X^{-E)\ X,iE) 



functions) . The spectral statistics of a given universality 
class can be described (and defined) by some ensemble 
of random matrices (usually there will be a lot of dif- 
ferent ensembles that share the same universal spectral 
statistics). In this paper we will focus on the ergodic 
universality classes that can be described by Gaussian 
ensembles of Hermitian matrices in each of the ten sym- 
metry classes. Note, that three chiral symmetry classes 
and the symmetry classes BD and Dili fall into various 
subclasses - as the universal spectral statistics is different 
in each of these subclasses they define different ergodic 
universality classes in the same symmetry class. As in 
the previous section we will only discuss one subclass in 
each of these cases. In the symmetry classes A, Al, All, 
C and CI there is one unique ergodic universality class. 
In the chiral classes Alll, DDI and CII we restrict to 
|p — ql =0 (see section III C 2|) . Finally, the classes BD 
and Dili have two subclasses (see section Hi C 3ll and we 
will restrict ourselves to the subclasses D (DD-even) and 
Dlll-even. 

Let us mention that apart from the ergodic universality 
classes there are a lot of other physically relevant univer- 
sality classes within each symmetry class. In random- 
matrix theory these correspond to ensembles which are 
not equivalent to the Gaussian ensembles. For instance 
ensembles of banded or sparse Hermitian matrices can 
describe quantum systems in a localized regime |l5j |. 

In Andreev systems more specialized random-matrix 
ensembles can describe the so-called hard gap in the 
quasi-particle excitation spectrum that appears when a 
small part of the boundary of a normalconducting chaotic 
billiard is coupled to a superconductor 16]. If no mag- 
netic field is applied the resulting combined electron-hole 
dynamics near the Fermi level is no longer chaotic and 
the system does not belong to an ergodic universality 
class. 



A. The fluctuating part of the density of states 

To reveal universality in the statistics of quantum spec- 
tra the system dependent mean density of states has to be 
separated. This is done by writing the density of states 
as a sum 

d{E) = ^5{E- Ei) = dwcy\{E) + 5d{E). (30) 



III. UNIVERSAL SPECTRAL STATISTICS 

In the previous chapter we have summarized the sym- 
metry classification of quantum systems. It is completely 
general. We have not yet related it to universal spectral 
properties. This will be done in this section. In each sym- 
metry class there are several universal regimes with re- 
spect to their spectral statistics. A universality class is a 
subset of a symmetry class which share the same spectral 
statistics (or at least some universal spectral correlation 



In presence of Kramers' degeneracy (symmetry classes 
All, CII, and Dili) we define the density of states such 
that every doubly degenerate energy is counted only once 
in the sum d{E) = 8{E — Ei). Let us also introduce 
a degeneracy factor 5, where g = 2 for systems with 
Kramers' degeneracy and else 5 = 1. 

In equation (|30|l the first part c?woyi(£') is the average 
density of states which may be obtained by counting all 
states Ei in an interval E— ^<Ei<E+^, then the 
number Nj of states in that interval divided by D/ is the 
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average density of states 

Nr 

dwcyliE) = (31) 

For this to be well-defined it is necessary to choose Ej 
self-consistently in range such that i. Nj ^ 1 which is 
equivalent to taking the energy interval much larger the 
the mean spacing Ej ^ AE = -r-^ — , and ii. Ej is small 
compared to compared to the scale on which the resulting 
rfwcyi changes. 

In systems that allow for a classical limit and one may 
consider the semiclassical regime. The scale Ej is then 
chosen classically small {Ej ^ as fi ^ 0) but large com- 
pared to the mean level spacing. Thus the average den- 
sity of states is well-defined in the semiclassical regime. 
It is given by Weyl's law 

dweyliE) = / ^^^^SiE ~ ffelass(p, q)) (32) 

where i7ciass(p, q) is the classical Hamilton function, and 
/ the number of freedoms. This equation shows that the 
average density of states defined by Weyl's law is system 
dependent and universal features can only arise due to 
the fluctuating part Sd{E). Note, that Weyl's law gives 
the mean density of states on scales much larger than the 
mean level spacing. In the presence of mirror symmetries 
the fluctuating part 5d{E) may contribute to universal 
features in the density of states on the scale of the mean 
level spacing. 

For classically chaotic (hyperbolic) systems the fluctu- 
ating part of the density of states is given by Gutzwiller's 
trace formula 17j as a sum over periodic orbits of the 
classical system 

M(i?)= V -^^„cos-^. (33) 

p.o. a 

Here, ta is the primitive period of the orbit (the time 
needed for a single traversal), Ac, = ° ^° ^ is the 

stability amplitude of the periodic orbit {M^'^ is the re- 
duced monodromy matrix and fi^ the Maslov index) and 
Wa = p dq is the (reduced) action. Note, that hyper- 
bolic chaos is a strong condition on a classical system - all 
periodic orbits are hyperbolically unstable and isolated 
in the energy shell. 

The energy scale for universal features is given by the 
mean level spacing AE = ^J; ^ . Introducing rescaled 
energies E — eAE one obtains a density of states 

d{e) = 1 + 5d{e) (34) 

for the unfolded spectrum. 

B. Gaussian ensembles of random- matrix theory 

Each ergodic universality class can be associated to a 
Gaussian ensemble of random matrices. Within one class 



the Gaussian ensembles differ only by the dimension of 
their matrices. The universal features of spectral statis- 
tics are extracted in the limit of large matrices. 

In each Gaussian ensemble the probability for a Hamil- 
tonian matrix 7i (with symmetries according to one of the 
ten symmetry classes) has the form 

P(H)dMH)-le-^'^«'dMH) (35) 

where A'^ is a normalization constant, A is an overall scale 
that fixes the mean level spacing, and the measure d/i(7i) 
is given by Y\ dRe TLij dim Tiy where the product runs 
over all independent elements of TL. 

In general, on may denote the Gaussian ensemble for 
the symmetry class X by X-GEi. We will use this notion 
for the Andreev classes. Note, that for some symmetry 
classes one should distinguish various ergodic universality 
classes. As we have restricted our investigations to just 
one relevant subclass we will use the name of the whole 
symmetry class for the Gaussian ensembles. 



C. Spectral form factors 

Let us now define the statistical functions that are in 
the center of our investigation. 

For a physical system the following averages are either 
performed over some system parameters or over different 
parts of the spectrum. We will always use unfolded spec- 
tra with unit mean level spacing. Spectral averaging is 
only possible if the universal results are invariant under 
shifts of the energy E ^ E + E'. 

There is an important difference between the Wigner- 
Dyson classes where the universality was conjectured for 
a single spectrum of one system and the remaining seven 
symmetry classes where some universal features near en- 
ergy E = can only be obtained by averaging over differ- 
ent spectra. In the classical limit one naturally obtains 
many spectra for the same physical system by formally 
changing h. Thus even for the seven novel symmetry 
classes one may average over different spectra for the 
same physical system. 

We will be interested in the two simplest correlation 
functions and their Fourier transforms. We will call the 
latter form factors. The first correlation function is sim- 
ply the averaged fiuctuating part of the density of states 
{Sd{e)). If the spectral statistics is invariant under shifts 
this expectation value must vanish (if not it would be 
a constant over scales much larger than the mean level 
spacing - in contradiction to its definition) . The spectral 
statistics near a spectral mirror symmetry is not invari- 
ant under energy-shifts. Non-trivial contributions to the 
mean fiuctuating part of the density of states may then 
arise. These have to appear on the scale of mean level 
spacing (else it would be inconsistent with the separation 
of the density of states in c?vvoyi + ^d) . 

The Fourier transform of the averaged fiuctuating part 
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of the density of states is the first- order form factor 



A'i(t) = 2 



{Sd{e)) 



(36) 



where tn = §^ is the Heisenberg time. Inverting the 
Fourier transform one may represent the deviations from 
Weyl's law in the expectation value for the density of 

states as 



(M(6)) 



/"OO 

/ dr cos (27rer) iiri(T). 



(37) 



Note, that for r > the first-order form factor is the 

trace of the time evolution operator Ki{t) — tre' s . 
The second-order correlation function is defined by 



C(e,eo) = {Sd{eo + e/2)Sd{eo ~ e/2)) 



(38) 



If the spectral statistics is invariant under energy shifts it 
only depends on the energy difference e - averaging over 
different parts of the spectrum for a given system is an 
average over eg. Its Fourier transform with respect to e 
is the second-order form factor 



dee 



'C(e,eo) 



(39) 



where we have suppressed the possible dependency on eg. 

For physical spectra a time average over a small time 
interval Ar <C 1 has to be added to the definition of the 
form factors. 



D. Spectral statistics for the Gaussian random 
matrix ensembles 



We will now summarize the relevant results from 
random-matrix theory (for more details see 0, 0, IE S 



I)- 



where Si(a;) = d^ is the sine integral. The corre- 
sponding second-order form factors are given by 




(42) 



|r| for |t| < 1 
1 for |t| > 1 

|r|(2-log(2|r|+l)) for |t| < 1 

2-k|logJ^ for|r|>l 

^(2-log||r|-l|) for|r|<2 

1 for |t| > 2 



2. The novel ensembles 

The Gaussian random-matrix ensembles in the chiral 
symmetry classes are known as chGUE (AIII-GE), ch- 
GOE (BDI-GE) and chGSE (CII-GE). The Andreev en- 
sembles C-GE, CTGE, ZJ-GE, and DIII-GE do not have 
any established name. The spectral statistics of these 
ensembles is not invariant under energy-shifts and, as a 
consequence, deviations from Weyl's law need not vanish 
near e = 0. At energies much larger than the mean level 
spacing |eo| ^ 1 Wigner-Dyson statistics is recovered. 
Thus, for the two-point correlation function we have 



(43) 



The universal features near the symmetry point e = are 
most prominent in the density of states. Though there 
are universal deviations from Wigner-Dyson statistics in 
all correlation functions we will focus on the density of 
states. The universal deviations from Weyl's law for the 
chiral ensembles are given by jj, ^ ilSj 11^ |23 



1. The Wigner-Dyson ensembles 

The ergodic universality classes for quantum systems 
in the Wigner-Dyson classes are described by the well- 
known Gaussian ensembles of random matrix theory 
GUE (^-GE), GOE (AI-GE) and GSE (^E-GE). The 
universal spectral statistics is invariant under shifts of 
the energy e — > e 4- eg- Thus the expectation value of the 
fluctuating part of the density of states vanishes and so 
does its Fourier transform 

Ki{t)^-^=Q. (40) 

The two-point correlation functions are given by 



^^ii(.;g(^.) + j2(,.))-i 

(5d'^"C'™(e)) + f Jo(-.)(l-Jo''' = '<^^.^(^)) 



(44) 



^(^) = /o 
:-M2 3.(^1.1)-.) (41) n(y, x) = /J 



The corresponding first-order form factors can be 
calculated explicitly in terms of the complete el- 
liptic integrals of first, second and third kind 

£{x) = x/T 



\J \ — x sin^ 



X sm 



and 



iSi(27r|e|) 



^ . ^ dcf) (we use the conven- 

(1— y sin^ a: sin^ 

tion that Il{y,x) is real for y > 1 [IJ). They are given 
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by 



41x1 



r( 41x1 

xl(l + k|) '^l, + 



21x| tt/ 4|t| 4|t| 

^TRTTTTaFTW 2F1+T ' -(T+TTip 



^(k|+2)(|r|+l)^ 



+ 1 ' (2+|r|)2 



+ 1 c-( 4|t| a |x|-1 4|t| \ 

W^l(|x| + 1)2 j+ '^l (Ixl + l)^;- 



(45) 



when an average over various spectra is performed. This 
corresponds to an average over some system parameter. 
We conjecture that for classically chaotic systems with a 
spectral symmetry all correlation functions of the fluc- 
tuating part of the unfolded density of states 6d{e) as 
given by (|34|l averaged over one system parameter coin- 
cide with those of the corresponding Gaussian random- 
matrix ensemble in the novel symmetry classes. This 
includes universal deviations from Weyl's law in the den- 
sity of states itself. Note, that though there are seven 
symmetry classes there are infinitely many ergodic uni- 
versal classes due to the different subclasses. Though 
some average is certainly necessary we may still conjec- 
ture the fidelity to ergodic random-matrix theories of a 
single physical system by formally averaging over spec- 
tra for different values of an effective Planck's constant. 
In superconducting-normalconducting hybrid structures 
this corresponds to an average over Fermi energy /i. 



= e(l-|r|) 



Hr) = iKi(5)Ci^GE^i_g(|^|_i)^ 



IV. QUANTUM STAR GRAPHS FOR THE TEN 
SYMMETRY CLASSES 



E. The generalization of the 
Bohigas-Giannoni-Schmit conjecture 

It has been conjectured by Bohigas, Giannoni and 
Schmit that quantum systems (in the semiclassical 
regime) with a chaotic classical limit have universal spec- 
tral fluctuations that coincide with the predictions of 
one of the Wigner-Dyson Gaussian ensembles of random- 
matrix theory GUE, GOE, or GSE. More precisely in an 
average over different parts of the unfolded spectrum the 
n-point correlation functions for n > 2 of a single spec- 
trum are conjectured to coincide with the correspond- 
ing correlation functions of the Wigner-Dyson ensemble. 
The mean density of states of a given quantum system is 
non-universal and cannot be described by random-matrix 
theory. Semiclassically it is given by Weyl's law. 

A lot of evidence has since been gathered both numer- 
ical and analytical that this conjecture is true in generic 
chaotic systems 5] (though a few exceptions are known 
[22,|2l|). Many approaches have been used to understand 
and proof the fidelity to random-matrix theory in com- 
plex quantum systems |3j HE HS] • Recently there has 
been considerable progress in the semiclassical approach 
using periodic orbit theory |27l |. 

Bohigas et al stated their conjecture before the impact 
of spectral mirror symmetries on spectral statistics has 
been recognized. A proper generalization of their state- 
ment has to take into account that a spectral average will 
wipe out all effects of a spectral mirror symmetry. Thus 
the original conjecture is expected to hold for the novel 
symmetry classes as well: after averaging over different 
parts of a single spectrum they will show the universal 
spectral fluctuations of GUE, GOE, or GSE. 

The additional universal features in physical systems 
near the spectral symmetry point can only be observed 



Quantum graphs have been introduced by Kottos and 
Smilansky |1Q|| as simple quantum systems with an exact 
semiclassical trace formula for the density of states. They 
consist of V vertices and connected by B bonds. Each 
bond bi connects two vertices and has a length Li. A 
particle propagates freely on the bonds and is scattered 
at the vertices by prescribed boundary conditions which 
leads to quantization. In their first approach Kottos and 
Smilansky considered vertex boundary conditions that 
implied current conservation and continuity of the wave 
function. The continuity condition is not always essential 
and has often been relaxed. In that case the boundary 
conditions at a vertex are specified by any unitary scat- 
tering matrix that transforms incoming waves to outgo- 
ing waves - unitarity of the vertex scattering matrix is 
equivalent to current conservation. 

We will not discuss general graphs but limit ourselves 
to a very simple class of graphs - star graphs. 



A. Quantization of star graphs 

A star graph consists of B bonds hj of length Lj [j = 
1, . . . , S) and V = B + 1 vertices vj {j — 0, . . . , B). Each 
bond bj emanates from the central vertex vq and connects 
it to the peripheral vertex Vj (see figure^. 

We will allow for a multi-component wave function on 
the graph. The number of components /i is assumed to 
be equal on all bonds. It may represent different spin 
components or electron and hole components of a quasi- 
particle. The /z-component wave function on the bond bj 
is 

*(^)(:r(^")) = e"=-"' + ct,lf> e'^^f^^— (46) 
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where < x^-') < Lj is the distance from the central 
vertex and 



^in (out) 



(47) 



are ^-component vectors of constant coefficients for the 
incoming (outgoing) waves on bond bj ("incoming" and 
"outgoing" wiU always be used with respect to the central 
vertex) . 

It is convenient to combine all the coefhcients of in- 
coming (outgoing) waves into two vectors of dimension 



4>u 



(out) 



.(B) 

''in (out), 1 



fin (out),M 



(48) 



The boundary condition at the center can then be written 
in the form 



'/'out =ScC(k) cj)i^. 
The diagonal /i_B x jiB matrix 



(49) 



(50) 



describes the propagation along the bonds. Here (and 
in the rest of this section) a = 1,2, .. ./i indicates the 
component of the wave function and j = 1, 2, . . . , i? is 
the bond index. The central scattering matrix Sc is a 
fixed unitary fiB x fj,B matrix that defines the boundary 
conditions at the center. For definiteness, we will assume 
that different components of the wave function do not 
mix at the center, thus 



c (a) iw\, • ., 



(51) 



where the B x B matrix S^'' describes scattering of the 
a-component at the center. 

The boundary conditions at the peripheral vertices 
may be described by one fixed x ^ vertex scattering 
matrices for each peripheral vertex - these can be com- 
bined to a single unitary fiB x fiB peripheral scattering 
matrix Sp such that 



SpC{k) cj)^ 



(52) 



Since different bonds are not coupled at the peripheral 
vertices 



5p,aj,a',' = = ^'apL'S'""^-"' (53) 



if) 




FIG. 1: Sketch of a star graph with five bonds of equal length. 



where cr(j) is the fi X fi scattering matrix at the vertex 

Vj. 

We will not allow any dependence of the scattering 
matrices Sc and Sp on the wave number k. Uniqueness 
of the wave function and the boundary conditions l|49l) 
and H52|l lead to the quantization condition 



</>i„ = SpC{k)ScC{k) (j)^^ = 5b (fc) 



(54) 



where we introduced the bond scattering matrix Spik). 
Non-trivial solutions of these equations exist only when 
the wave number belongs to the discrete spectrum k — kn 
given by the zeros of the corresponding determinant 

det(5B(fc„) - 1) = 0. (55) 

The density of states for the graph is defined as 



oo 

d(A:) = - V<5(fc-fc„), 

^ n=0 



(56) 



where 5 = 2 in systems with Kramers' degeneracy (else 
a = 1). 



B. The trace formula 

Let us now write the density of states as a sum of its 
mean dwcyi and an oscillating part 6d{k) 



d{k) = dwoyi + Sd{k). 



(57) 



For both contributions one can give an exact semiclassical 
expression. The mean density of states is given by Weyl's 
law 



iWoyl 



gn 



(58) 



and the oscillating part obeys the trace formula ^3 

7 1 

5d(k)^lm — Y tr5s(fc)". (59) 



dk ^ — ' qiT n 



In the sequel we will consider star graphs where all bond 
lengths are equal Lj — L. In that case the bond scatter- 
ing matrix is a periodic function of k 



SB{k)=SBik + j-) = e'^'''^SB 



(60) 
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where 



SB^SB{k = 0)=SpSc 



(61) 



is the reduced bond scattering matrix. Thus, the spectrum 
is also periodic and the trace formula simplifies to 



d{k) = 



fiB L 2L 



gn 



giT 



n=l 



(62) 



The periodicity of the spectrum will not be relevant here 
as we are interested in features on the scale of a mean 
level spacing. 

For equal bond lengths the trace formula can be de- 
rived in a few lines: Let e~"^^ [j = l,...,fj,B) be the 
eigenvalues of the unitary reduced bond scattering ma- 
trix Sb ■ The quantization condition (|55|l is equivalent to 

^ ~ li n^od r (j = 1: ■ ■ • J M-S) and the density of states 
is 



j—1 71 — — OQ 



3 = 1 



(63) 



n— — OQ 



The mean density of states d-yvcyi = ^^^^ is just the 
n = term in the sum over n while the rest gives the 
trace formula for the oscillating part (the second line 
follows from the first by Poisson's summation formula 



The trace formula (|59|l can be interpreted as a sum 
over periodic orbits p on the graph. A periodic orbit p = 
[(ji, ai), (j2, a2), • • ■ (jn, Q;„)] of length n is defined by a 
sequence of n peripheral vertices Vj-^ Vj^ . . . Vj^ visited one 
after the other together with the specification of the wave 
component aj between two vertices (cyclic permutations 
define the same orbit). A periodic orbit is primitive if it 
is not the repetition of a shorter periodic orbit. In terms 
of primitive periodic orbits p and its repetitions the trace 
formula reads 



expression for the form factors that can easily be aver- 
aged numerically. In the second paper of this series the 
trace formula will be in the center of an analytic approach 
to universality. 

Since universality exists on the scale of the mean level 
spacing we will write k — AtAfc where Afc — -j^^ is the 
mean level spacing. In terms of the rescaled wave number 
the trace formula is 



d{K) = l + Re y e'^^'^^s.^. (65) 

' n=l 

We have introduced the shorthand 



= tr53 



(66) 



for the 7i-th trace of the reduced bond scattering matrix. 

The first-order form factor is obtained by a Fourier 
transform and a subsequent time average. It obeys the 
trace formula 



Ki(t) 
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gn 

]Ib 



Ki, 



(67) 



where the bar denotes a time average over a small time 
interval At = ^ < 1 and 



g 



(68) 



The brackets (•) denote an average over an ensemble of 
graphs. This can be written more compactly as 



Ki{t 



fiB' 



Ki. 



1 
An 



An-l 

E 

k=Q 



Ki,n+k 



(69) 



where the continuous time average has been replaced by 
an average over the discrete time t = 

The second-order form factor for a graph also obeys 
a trace formula which, after a spectral average over the 
central wave number, is given by 



2L 

5d{k) = — np (Api 

^ p.p.o.:p r— 1 



(64) 



where 



is the length of the primitive periodic 



orbit — 



apai + iaiClc 



Jl + lJl 



is the ampli- 



tude of the primitive orbit and Wp 



ji+iji 



~r aipoii^iai 



its phase ("action"). Note, 

that we set jn^+i = ji and a„p+i = ai. 

The similarity of the sum over periodic orbits 164() 
to the semiclassical Gutzwiller trace formula is evident. 
However, while semiclassics, in general, is an approxima- 
tion the semiclassical trace formula for quantum graphs 
is exact. 

The trace formula will be our main tool in the analysis 
of universal spectral statistics. It will lead us to a simple 



where 



K2,a = 



gfiB 



(70) 



(71) 



If no spectral average is performed additional terms ap- 
pear. These are irrelevant for the graphs in the Wigner- 
Dyson classes (they do not survive the subsequent ensem- 
ble average) . Here, we will not consider the second-order 
form factor for graphs in the novel symmetry classes 
where the additional terms are relevant near the central 
energy eo = 0. 

Though Ki^n and K2,n do not involve a time average 
we will refer to them as (discrete time) form factors. 
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C. Star graphs for all symmetry classes 

We will now construct ensembles of star graphs for each 
symmetry class. The star graphs will be constructed in 
such a way that spectral fluctuations of the corresponding 
ergodic universality classes can be expected. Though we 
are not able to prove an equivalent conjecture we will 
give strong evidence. 

The constructions of star graphs for each symmetry 
class are based on a proper choice of the central and 
peripheral scattering matrices Sc and Sp. Both have to 
obey the right symmetry conditions (see section Hljl . 

As an example of star graphs in class Al that do not be- 
long to the ergodic universality class let us mention Neu- 
mann star graphs. These have a one-component wave 
function (/i = 1) on the bonds, Dirichlet boundary con- 
ditions at the peripheral vertices such that Sp = — 1 
and Neumann boundary conditions at the center, thus 
Sc,ki =_s.~ ^ki- Such graphs have been investigated 
first in 13 and in more detail in - in contrast to 
our approach the bond lengths were chosen different for 
each bond (and incommensurate). However, Neumann 
boundary conditions at the center favor backscattering 
and lead to non-universal (localization) effects psf. 

Our approach is different in as much as we allow for 
more general scattering matrices at the center and in as 
much as we will always consider an ensemble of graphs. 
The occurrence or non-occurrence of localization effects 
can be traced back to a gap condition on the matrix 
Tij = \Sc-ij\'^- This bistochastic matrix describes the 
corresponding "classical" probabilistic dynamics (equiv- 
alent to a Frobenius-Perron operator). In chaotic (er- 
godic) systems the Frobenius-Perron operator has a finite 
gap in the spectrum between the (unique) eigenvalue one 
and all other eigenvalues which describe the decay of the 
probability distribution. In Neumann star graphs this 
gap is small and vanishes in the limit i? — > cxd faster than 

which leads to non-ergodic spectral statistics. 

In general one needs a multi-component wave function 
to introduce the different symmetries. The number of 
components /i has been chosen minimal under the addi- 
tional assumptions that the components do not mix at 
the central vertex and that time-reversal is only broken 
at the peripheral vertices. 

Though we explicitly choose the central and peripheral 
scattering matrices guided by simplicity and minimality, 
most of the results are much more general. 

The central scattering can be chosen in a very sim- 
ple way by using the symmetric B x B discrete Fourier 
transform matrix (29l] 

SuFT.ki = 4=e''"^ (72) 
V B 

or its complex conjugate for each component. An incom- 
ing wave on a given bond is scattered with equal prob- 
ability to any bond which excludes localization effects. 
Indeed, the matrix TijFT.fci = |'5DFT,fciP = -g has one 



eigenvalue 1 while all other B — 1 eigenvalues vanish (so 
the gap is maximal). 

The bond scattering matrix 5b (fc) = ScC{k)SpC{k) 
for each ensemble of graphs is constructed by demand- 
ing that the matrices Sc, Sp, and C{k) do all have the 
canonical forms of the desired symmetry class given in 
section m Note, that for star graphs we are interested 
in the k spectrum, so in the canonical forms for scatter- 
ing matrices the energy E has to be replaced by k. The 
ensemble of graphs is built by introducing some random 
phases into the peripheral scattering matrix. 

1. Star graphs in the Wigner- Dyson classes 

Let us start with the simplest case: an ensemble of star 
graphs in class Al where a one-component wave-function 
suffices to incorporate the time-reversal symmetry which 
demands that the unitary matrices Sc, Sp, and C are all 
symmetric. Now, Sp and C{k) = e''^'^l are diagonal for 
/i = 1, and choosing 

Al: Sc^ 5dft (73) 

we meet all requirements. At the peripheral vertices we 
are free to choose one random phase Pk for each periph- 
eral vertex j independently such that 

Al : Sp^ki = Skie"^' (74) 

where < /3fe < Stt is uniformly distributed. 

For class A we have to break time-reversal symmetry. 
This may be done by choosing a non-symmetric central 
scattering matrix for a one-component wave-function on 
the graph. As we like to keep the simplicity of the dis- 
crete Fourier transform matrix we choose another simple 
construction with a two-component wave function, a cen- 
tral 2B X 2B scattering matrix 

A: Sc= f ° ] , (75) 
y 5dft J 

and 

A: Sp = ±h ""A (76) 

for the peripheral scattering matrix. The diagonal ma- 
trices Vj are 

Vi,M = (5fc,e'('3^+T'=) 

where the independent random phases f3k, 7fe, and Sk are 
uniformly distributed. 

For class All a four-component wave- function is needed 
to incorporate time-reversal invariance with = -1 
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FIG. 3: First-order form factor for Andreev star graphs in 
the classes C and CI averaged over 10000 realizations with 
B = 100 bonds (additional time average over an interval of 
length 8-^). Dashed lines: prediction by Gaussian random- 
matrix tfieory. Full lines: numerically calculated form factor 
for graphs. Inlets: discrete time form factor as function 

of r = n-^ for even n {Ki,n vanishes by construction for odd 



FIG. 2: Second-order form factor for star graphs in the 
Wigner-Dyson ensembles averaged over 10000 realizations of 
with B = 100 bonds (additional time average over an in- 
terval of length 8^). Dashed lines: prediction by Gaus- 
sian random-matrix theory. Full lines: numerically calculated 
form factor for graphs. Inlets: discrete time form factor K2,n 
as function of r 
the class All graphs 

2 



g for the class A and Al graphs. For 
2,n vanishes for odd n - the inlet shows 



for even n. 



into our scheme for star graphs. Indeed, the number of 
components must be even as discussed above in III C II 
In addition, we assumed that components do only mix at 
the peripheral vertices. Then, a 4 x 4 scattering matrix 
at the peripheral vertices is the minimal matrix dimen- 
sion that allows for component mixing as can be seen 
from the canonical form H16() of a All scattering matrix. 
The diagonal matrix C{k) is a diagonal unitary matrix 
of the canonical form. All further requirements are met 
by choosing 



All 



Sc 



/5dft \ 

5dft 

5dft 

V 5dft/ 



(78) 



for the central 4B x AB scattering matrix, and 



All : 



1 

71 



f Vi 
V3 

-V4 
\Da 



V2\ 

~V2 

I?3 

Vi / 



(79) 



for the peripheral scattering matrix. The diagonal ma- 
trices I?j are given by 



All : 



Vim = 4ie'('3'=+'^'=) 
■Dam ^ dkie'^f^-'-l 



(80) 



where the independent random phases f3k, 7fc , and Sk are 
uniformly distributed. 

Ergodic spectral statistics may be expected for the 
three ensembles in The Wigner-Dyson classes. This is 
strongly supported by a numerical calculation of the 
second-order form factor (see figure O. 



2. Chiral and Andreev star graphs 

Let us start with the Andreev star graphs for the 
classes C and CI where the wave function can be cho- 
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FIG. 4: The first-order form factor for Andreev star graphs 
in the classes D and Dili (see figure for details). 



sen in the simplest case to have two components. The 
first will be called "electron" and the second "hole" . The 
transfer matrix C{k) and the central scattering matrix 
defined by 



Sc 







5, 



(81) 



DFTy 



obey the symmetry condition (|27|l . 

The peripheral scattering matrix may be chosen such 
that complete Andreev scattering (electron-hole conver- 
sion) takes place 




Sp = 

where the diagonal matrix V is 
for class C, and 

CI : T>ki = Ski io-fe 



(82) 



(83) 



(84) 



for class CI. The random phases (3k are uniformly dis- 
tributed and (jfc = ±1 with equal probability. 

For the Andreev classes D and Dill and as well for the 
three chiral classes Alll, BDl, and CII a four-component 
wave function is needed. We will call the first (last) two 
components "electron" ("hole"). 

The symmetry requirements , (1^ , Ip^ , |2U , and 
are met by the transfer matrix C{k) and the central 



ifl(r) 



class ,4111 (chGUE) 




Kl{r) 



class BDl (chGOE) 



Klir) 




FIG. 5: The first-order form factor for chiral star graphs (see 
figure |3] for details). 



scattering matrix defined by 

D 

Dili 

„ U i^DFT U U . ^ 

All! > : oc = „ • (85) 

I 5^,^^, ' 

BDl 
CII 

In all five remaining classes we choose the peripheral 
scattering matrix such that complete Andreev scattering 
takes place. For D and Dill the simplest choice obeying 
the symmetry requirements are 











\ 





5dft 














'^DFT 





\ 








'^DFT/ 




1 

71 



( 


V 



Vi V2\ 
-V^ VI 
V* 

0/ 



1 -^2 

\-V2 Vi 



(86) 



where the diagonal matrices "Dj are 
T), , , 

D : 



1^2, kl 



Ski e'^" 
Ski e'^" 



(87) 
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for class D, and 



Dili 



1^2. kl 



Ski e'^' 

Ski iCTk 



for class Z5III. The random phases (3k and are uni- 
formly distributed and crfc = ±1 with equal probability. 

The simplest choice for peripheral scattering matrices 
in the chiral classes is 



(89) 



The diagonal matrices Vj have to be chosen according to 
the requirements of each symmetry class. For class Alll 
they are 






/o 







V2 \ 


1 














Vi 














-I?4 





/ 



Alll : 



■Dim 


= SkiCTk 




= Skie'"" 




= Skie-'^" 




= SklTk 


2^5, fei 


= Skie'^" 


2^6, fci 


= Skie-'^" 



(90) 



where Tk,ak = ±1 with equal probability and the phases 
(3k and are uniformly distributed. The BDl star 



graphs can be obtained from the AIII case by the ad- 
ditional restrictions 



BDl 



Tk = CFk 



and 



Ik 



-(3k- 



(91) 



Finally, for class CII the peripheral scattering matrix is 
defined by 



CII : 



= SklUk 

= Ski&f^" 

= —SklCTk 

= -Skic'^'- 
- -Skie-""" 



(92) 



We have checked numerically that the first-order form 
factor for the constructed Andreev and chiral star graphs 
obeys the corresponding prediction by Gaussian random- 
matrix theory (see figures I^Bl andjSJ- 
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